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^ ■ ABSTRACT 

The Spitzer spectrum of the giant FR II radio galaxy 3C 326 is dominated 

^ ■ by very strong molecular hydrogen emission lines on a faint IR continuum. The 

^ ' H2 emission originates in the northern component of a double-galaxy system 

Q ■ associated with 3C 326. The integrated luminosity in H2 pure-rotational lines is 

■ 8.0 X 10^^ erg s~^, which corresponds to 17% of the 8-70 /xm luminosity of the 
^ ■ galaxy. A wide range of temperatures (125-1000 K) is measured from the II2 0-0 
O ■ S(0)-S(7) transitions, leading to a warm H2 mass of 1.1 x IO^Mq. Low-excitation 

ionic forbidden emission lines are consistent with an optical LINER classification 
^ ■ for the active nucleus, which is not luminous enough to power the observed H2 

■ emission. The H2 could be shock-heated by the radio jets, but there is no direct 
indication of this. More likely, the H2 is shock-heated in a tidal accretion flow 
induced by interaction with the southern companion galaxy. The latter scenario is 
supported by an irregular morphology, tidal bridge, and possible tidal tail imaged 
with IRAC at 3-9 /im. Unlike ULIRGs, which in some cases exhibit H2 line 
luminosities of comparable strength, 3C 326 shows little star-formation activity 
(~ 0.1 M0 yr~^). This may represent an important stage in galaxy evolution. 
Starburst activity and efficient accretion onto the central supermassive black hole 
may be delayed until the shock-heated H2 can kinematically settle and cool. 



Subject headings: galaxies: active, galaxies: jets, infrared: galaxies 
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Introduction 



Recent Spitzer and ISO observations liave revealed a new class of extremely luminous 
molecular hydrogen emission galaxies (MOHEG's) with L{B.2) = 10^° — 10^^ erg s~^ in 
pure-rotational molecular hydrogen emission lines but relatively weak total IR emission, 
L(H2)/L(IR)> 10^^. In comparison, normal sta r-forming galaxies in the SINGS survey of 
the local universe have L(H2)/L(IR)~ 6 x 10"^ JRousscI et al~ll2007h . 



Applet on et al. I (120061 ) discovered H2 emission extended over 50" in the giant intergalac- 
tic X-ray shock front in the Stephan's Quintet compact group. The estimated H2 luminosity 
from the whole shock is in excess of 10^^ ergs/s. Nearly 30% of the IR lumin osity is in H2 
lines, while the rest is emitted by ~ 40 K dust in the far-IR (jXu et al. II2003I ). The power 
in the H2 lines compared with the dust continuum is several orders of magnitude larger 
than that seen in photo-dissociation regions (PDRs). The H2 emission lines are very broad 
(870 km s~^), suggesting turbulent motions inside the shock front. The H2 emission may be 
powered by non- dissociative shocks in dense clumps or filaments inside the shock front. 



Egami et al. I (120061 ) find extremely luminous L{B.2)= 2 x 10^^ erg s ^ emission from 



~ 10 M0 of warm H2 in the brightest galaxy of the high-luminosity X-ray cluster Zw 3146. 
This is accompanied by high-luminosity dust emission powered by starburst activity. Lower 
l uminosity ro- vibrational H2 emission is similarly seen in several l ow-redshift galaxy clusters 
( iDonahue et al. II2OOOI : iHatch et al. II2005I : I Johnstone et al. 1120071 ). The source of excited H2 
in galaxy clusters is an unsolved mystery-cooling flows, jet shocks, and X-ray heating by 
active galactic nuclei (AGNs) have all been considered. 

Molecular hydrogen emission is often seen in interacting galaxies and ultra-luminous 
infrared galaxies (U LIRGs), but L(H2)/L(IR) is typically « lO'^ Jflaas et al. 1 12005 



Higdon et al. II2006I ). The late-stage merger/LIRG NGC 6240 has atypically luminous emis- 



(Lutz et al. 


V\ CLl 111 

2003; 


j.j.2i ai^i^uiiipcii 

Netzer et al. 


2005; 


Max et al. 


2005; 


Armus et al. 


2006) 



The overlap region between the Antennae colliding galaxies also shows unusual ly strong H2 
emiss ion with L(H2)/L(IR)~ 1 x 10~^, possibly indicating a pre-starburst shock (IHaas et al 
2005h . 



We present the remarkable H2 emission spectrum of the FR II radio galaxy 3C 326, 
observed with the Spitzer Space Telescope. This is part of a large sur vey of 52 FR II radio 
galaxies from the 3CRR catalog (lOgle. Whvsong. fc Antonucci 1 120061 ) . and another survey 
of 21 FR I radio galaxies ( Ogle. Antonucci. fc Whysong I 2007 ). So far, we have found 3 FR 
lis (6%) and 5 FR Is (24%) with strong H2 emission (Ogle et al. 2007, in preparation). Of 
the galaxies in our sample, 3C 326 has the largest H2 emission line equivalent widths. The 
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weak mid-IR continuum of 3C 326 indicates the presence of a radiatively inefficient AGN, in 
spite of its large radio power and enormous size of its radio lobes. 



2. Radio Galaxy Identification 



The radio source 3C 326 is one of the largest known (jWillis fc Strom I Il978l ). with 
an angular size of 19'.5 (1.9 Mpc; we assume a cosmology with Hq = 70 km s~^ Mpc~^, 
Qrn = 0.3, and = 0.7). The pair of galaxies 3C 326 N and 3C 326 S are both at a 
redshift of z = 0.089 ±0.001 and have a projected separation of 24'.'9 (41 kpc). There is some 
confusion concerning which of the two galaxies hosts the AGN that produced the Mpc-scale 
radio lobes. Both galaxies have unresolved (< 3'^ ) radio sources coincid ent with their nuclei, 



with S'y(8 GHz) = 1.5 and 6.5 mJy, respectively (IRawlings et al. Ill990l ). While the southern 
radio source is brighter, both sources are consistent with typical FR II radio galaxy cores. 
The galaxy 3C 326 N has a LINER-like optical emission line spectru m (strong [On], Ha, [S 



III, and weak [C 


nil), 


Simpson et al. 


1996) 



Also, the K-band magnitude and corresponding ste llar mass of 3C 



326 N is more typical of FR II radio galaxy hosts (IRawlings et al. 
targeted 3C 326 N with the Spitzer IRS. 



1990). We therefore 



3. Observations 
3.1. Spitzer IRS 



The radio galaxy 3C 326 N was observed with the Spitzer IRS ( iHouck et al. II2004J ). 
covering the wavelength range 5-35 /xm. The companion galaxy 3C 326 S was observed on 
one nod of the Short-Low (SL) slits because of a fortuitous slit PA (Fig.l). The slit widths 
in the low-resolution modules Short-Low 2 (SL2), Short-Low 1 (SLl), Long-Low 2 (LL2), 
and Long-Low 1 (LLl) are 3'.'6, 3'.'7, 10'.'5, and 10'.'7, respectively. The SL and LL spectral 
resolutions range from 60-130. The exposure times were 240 seconds in each spectral order 
except LL2 (480 seconds). 

Our data reductions began from the S15.3.0 pipeline-processed basic calibrated data 
(BCDs). Off-slit observations were subtracted to remove background light. The 2D-spectra 
were cleaned using IRSCLEAN 1.8 to median filter bad pixels. Spectra were extracted using 
SPICE 1.4 within tapered regions matching the 5'^}2teer point-spread function (SL2: 7'.'2 at 
6 /im, SLl: 14'.'4 at 12 jjm, LL2: 21'.'7 at 16 jjm, LLl: 36'.'6 at 27 jjm). The second LL2 nod 
was considerably noisier than the first, so we discarded this data. Optimal extraction (which 
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assumes a point-source spatial profile) yielded significant gains in S/N for both galaxies. 

At the redshift of 3C 326, I'.'O corresponds to 1.7 kpc. The SL and LL slit widths 
correspond to 6.2 kpc and 18 kpc, respectively, and are considerably smaller than the optical 
extents of the two galaxies. However, the continuum level is well-matched across orders and 
the emission lines and dust features do not appear to be spatially extended along the slit 
for either galaxy (Figs. 1, 2, 3). The continuum and emission line fluxes from the optimally 
extracted spectra also agree well with the regular extractions, indicating that the galaxies 
are (excepting the stellar component at short wavelengths) mostly unresolved by Spitzer 
IRS, and that aperture effects are unimportant for studying the dust and gas emission. 



3.2. The Mid-IR Spectra 

The spectrum of 3C 326 N is dominated by pure-rotational emission lines from molecular 
hydrogen (Fig. 2, Table 1). We detect the full series of emission lines from H2 0-0 S(0)-S(7). 
Forbidden emission lines from the low-moderate ionization species Fe ll, S IV, Ne ll, S ill, 
and O IV are also observed (Table 2). Polycyclic aromatic hydrocarbon (PAH) emission 
features are seen at 6.7, 7.7, 10.7, 11.3, 13.6, and 17 /im (Table 3). We consider the PAHs 
at 6.7, 7.7, 10.7, 15.9 and 18.9 /xm and the Fe ll, S IV, and S ill lines to be marginal (< 3a) 
detections. 



We model the 3C 326 N spectrum using the fitting routine PAHFIT (jSmith et al. 



20071 ). Extinction is not included in the model since it does not improve the fit. Emission 
lines are fit by Gaussians and PAH features are fit by Drude profiles (Tables 1-3). None 
of the emission lines are resolved, all with FWHM< 2300 — 5000 km s^^. The 4.8-8 fim 
continuum is dominated by galaxian starlight. An excess over the stellar component in the 
8-24 fim region is fit by thermal emission from warm (90-300 K) dust, with integrated flux 



Fs-24 = 9.0 ± 0.9 X 10^^^ erg s"^ cm'^ and luminosity L8„24 = 1-8 ± 0.2 x 10^^ 

The companion galaxy 3C 326 S has no H2 emission detected; the 3a upper limit to the 
H2 S(3) line flux is 1.2 x 10"^^ erg s"^ cm^^ (l(H2)< 2.4 x 10^° erg s"^). No ionic forbidden 
emission lines are detected, consistent with the lack of optical emission lines and indicating 
an inactive nucleus. A spike in the spectrum at ~ 12 fim does not match any known emission 
lines and may result from a hot pixel. The 7.7 and 11.3 PAH features are 2-3 times as strong 
in 3C 326 S compared to 3C 326 N (Table 3) and the mean 9-10 /im continuum flux density 
(0.62 ± 0.04 mJy) is 20% brighter. 
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3.3. Spitzer IRAC and MIPS 



We retrieved the IRAC flFazio et al. Il2004f ) 3.6, 4.5, 5.8, and 8.0 fim S14.0.0-processed 



images from the Spitzer archive (Program-ID 03418, PI M. Birkinshaw) to further constrain 
the spatial distribution of the stellar, dust, and molecular gas components (Figs. 3,4). A 3a 
clipped average of the dither positions, excluding frames affected by a bright star at the edge 
of the field, was used to subtract the background light, and frames were combined using the 
SSC MOPEX-MOSAIC software. The 3.6 and 4.5 /im bands measure the stellar continuum. 
The 5.8 /im band contains the redshifted H2 S(5)-S(7) lines and 6.7 /im PAH feature. The 
8.0 /im band contains the redshifted H2 S(3)-S(4) lines and 7.7 /im PAH feature. Both of 
the latter two bands also contain thermal continuum emission from warm dust. 

The irregular morphology of 3C 326 N in the IRAC images suggests an ongoing inter- 
action or merger with 3C 326 S (24'.'9 at PA= 177°), a relatively compact elliptical galaxy. 
The isophotes of the two galaxies merge and they appear to be connected by a tidal bridge. 
There is a possible tidal tail extending from the NW quadrant at PA= 290 — 334°. Two 
other sources, one 4'.'3 away from the center of 3C 326 at PA= 111° and one 15'.'9 away at 
PA= 268°, appear to be unrelated background galaxies seen in projection, based on their 
morphologies, red colors, and small angular diameters in archival Sloan Digital Sky Survey 
(SDSS) images 0. 



We also obtained 24, 70, and 160 /im MIPS ( iRieke et al. 1 12004| ) imaging data from 
the Spitzer archive for 3C 326 (Fig. 5). Individual BCD frames from the S14.4 pipeline 
were processed with the MOPEX-MOSAIC software to generate final images. Aperture 
photometry was performed on both 3C 326 N and 3C 326 S, and point-source aperture 
correctionqj were applied to derive the fluxes (Table 4 and Fig. 6). Both galaxies are 
detected at 24 and 70, but not 160 /im. We measure a root-mean-square noise of 0.4 MJy 
sr~^ and 3a flux upper limits of 34 mJy at 160 /im. This is consistent with a non-detection 



at 870 /im (jOuillen et al. II2003I ). within a 23" beam centered on 3C 326 S. 



The companion 3C 326 S is twice as bright as 3C 326 N at 24 /im, while 3C 326 N is 
40% brighter at 70 /im. A third source is seen in the 70 /im image at the location where 
the tidal bridge in the IRAC 3.6 /im image connects to 3C 326 N (Fig.5, Table 4). It is not 
detected in the MIPS 24 /im band, nor in any of the IRAC bands. If it is truly associated 
with 3C 326, then it may indicate cool dust in the tidal bridge. Unfortunately, none of the 
IRS slits covered this feature, so we do not know if there is any associated H2 emission. 



^http://www. sdss.org/ 

^obtained from http://ssc.spitzer.caltech.edu 
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The spectral energy distributions (SEDs) of 3C 326 N and S are similar, but do show 
interesting differences (Fig. 6). Both are dominated by stellar emission at frequencies > 
4 X 10^^ Hz (wavelengths < 8 /xm), consistent with our model of the IRS spectrum of 3C 326 
N. Dust emission contributes most of the fiux at ~ 1 x 10^^ — 4 x 10^^ Hz and synchrotron 
emission from the radio cores is seen at 3-10 GHz. The ratio L(70 /im)/L(24 /xm) is greater 
in 3C 326 N than in 3C 326 S (Table 4), indicating relatively more cool dust emission and 
relatively less warm dust emission. 

We estimate far-lR continuum luminosities by interpolating the 24 and 70 fim photomet- 
ric points with a power law. This yields L24-70 = 3.0 ±0.4 x 10^^ erg s^^ and 3.3 ±1.8 x 10^^ 
erg s~^, respectively, for 3C 326 N and 3C 326 S. Adding together the MIR 8-24 /xm and FIR 
24-70 urn emission, the 8-70 pim luminosity of 3C 326 N is Lg-ro = 4.8 ± 0.5 x 10^^ erg s^^ 
(1.2 ± 0.1 x IO^Lq). Upper and lower bounds to the FIR luminosity at longer wavelengths 
of 

-^70-1000 — 0.3 — 1.7 X 10^^ erg s ^ are estimated for either galaxy by assuming blackbody 
emission that peaks at either 70 //m or 160 //m (Fig. 6). The total 8-1000 /xm IR himinosity 
of 3C 326 N therefore hes in the range Ljr = 0.8 - 2.2 x 10^^ erg s"^ (2 - 6 x IO^Lq). 



4. H2 Emission 

The galaxy 3C 326 N has H2 pure-rotational lines of extraordinary luminosity and 
equivalent width (Fig. 2, Table 1). The H2 S(3) hue has a luminosity of 2.34 ± 0.08 x 10^^ 
erg s~^ and an equivalent width of 0.80 ± 0.08 /im. The integrated S(0)-S(7) line luminosity 
is L(H2) = 8.0 ±0.4 X 10"^^ erg s^^. This is 17 ±2% of the measured 8-70 fim. luminosity and 
4 — 10% of the estimated 8-1000 fim luminosity, making it the most extreme H2-emitting 
galaxy seen so far by Spitzer (Table 6). Unlike the giant, extended, intergalactic shock 
in Stephan's Quintet, the H2 emission here is unresolved and contained within the central 
region of 3C 326 N. 

We estimate the upper-level column densities (N/g) for each observed H2 transition, 
divided by the statistical weights (Fig. 7, Table 1), assuming that the H2 source just fits 
inside the SL slit (3'.'7). The column density distribution is fit by a minimum of 3 temperature 
components with T — 125, 400, and 1000 K (Table 5). The fit is not unique but rather serves 
to demonstrate the large range of temperatures of the H2 emission regions. The total column 
density of 125 K H2 averaged over the SL slit is 1.8 x 10^^ cm~^. 

The ratio of ortho (odd angular momentum quantum number J) to para (even J) H2 
depends on the gas temperature and thermal history. Warm H2 in thermal equilibrium at a 
temperature of 125 K should have ortho/para=2.1, which is consistent with our model fit. 
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The ratio of ortho/para H2 in the hot (400 K and 1000 K) compone nts is consistent with the 



standard value of 3. for gas in thermal equilibrium at T > 250 K (jWilgenbus et al. 112000 



Neufeld et al. Il2006[ ) 



The masses of warm and hot molecular gas sum to 1.1 x 10^M(7^ (Table 6). This is 
comparable to the warm H2 mass in NGC 6240 (ILutz et al. II2003I : lArmus et al. II2006I ) and to 
the tota/ molecular gas mass of the Milky Way (~ 2 x lO^M©). No published CO observations 
of 3C 326 are available, so we have no knowledge of the mass of cold H2; although by analogy 
with galaxies like NGC 6240 the cold component could be even larger than the warm and 
hot components. In thermodynamic equilibrium, the H2 co mponent densities must exceed 
the critical densities of the observed rotational transitions (ILe Bourlot et al. Ill999l ). This 
yields maximum filling factors ranging from /^ax = 2 x 10~^ for the 125 K component to 
/max = 1 X 10^^ for the 1000 K component within the central 6.2 kpc-diameter sphere (Table 
5). 



5. Forbidden Emission Lines and AGN Activity 



The mid-IR forbidden emission lines from several ions (Table 2) can be used to distin- 
guish between starburst and AGN activity. In particular, line flux ratios of [O lv]/[Ne ll] 



0.6 ± 0.1 and [S IV 



( ISpinoglio fc Malkan 



/[Ne 11] = 0.08 ± 0.04 in 3C 326 are characteristic of a LINER AGN 



19921 : ISturm et al. II2006I ). Seyfert galaxies have relatively stronger [O 



iv] and [S iv] than this, while starburst galaxies are relatively weaker in these two lines. 
Consequently, little of the forbidden line emission comes from star-formation activity. A 
LINER classifica tion for 3C 326 is also consistent with the optical line ratio [O ill] 5007A/[O 



II] 3727A= 0.30 fISimpson et al. Ill996h . 



S'pzfeer observations give an upper limit to the IR luminosity of the AGN, since there is a 
significant contribution from the host galaxy. The fiducial 15 / xm continuum flux is 0.87±0.09 



mjy, somewhat greater than the upper limit measured by lOgle. Whysong. &: Antonucci 
( 120061 ). (Our measurement accuracy has improved by throwing out the noisy LL2 nod and 
using optimal extraction.) The corresponding luminosity is i/L(15 /im) = 3.4 ± 0.4 x 10^^ 
erg s~^ and its ratio to the radio lobe luminosity is only z/L(15 ;um)/z/L(178 MHz)= 4.4. 
We confirm that this powerful radio galaxy contains an AGN with relatively low accretion 
power. In comparison, we find from archival Spitzer IRS observations that the LINER/FR I 
radio galaxy M 87 is one tenth as bright, with z/L(15 fim) = 3.3 ±0.1 x 10^^ erg s~^, and has 
a similar ratio of z/L(15/im)/i/L(178MHz) = 3.9. It has been argue d that low-power radi o 



galaxies such as M 87 are fueled by radiatively inefficient accretion (IReynolds et al. Ill996l ) 
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X-ray he ating has been proposed as a possible mechan i sm to produce Hg emissi on in ac- 
tive galaxies (IKrolik fc Lepp lll989l : lDraine k. Woods lll990l : lRigopoulou et al. 1120021) . Unfor- 



tunate ly, there are no published X-ray observations of 3C 326 N. However, iHardcastle. Evans, fc Croston 
( 120061 ) find a ~1:1 correlation between 15 yum continuum emission and accretion-powered 2- 
10 keV X-ray luminosity for z < 0.5 radio galaxies. If this correlation holds true for 3C 326 N 
and the H2 were powered by X-rays, then the observed ratio of L(B.2) / L{15 /im)= 0.23±0.03 
would require an implausibly high conversion of X-rays to H2 emission. Chandra X-ray ob- 
servations are necessary to directly confirm this, and will also be useful for measuring the 
high-energy SED of the AGN to test accretion models. 



Aromatic Feature Emission and Star Formation 



We detect several aromatic (PAH) emission features in the IRS spectrum of 3C 326 
N (Fig.2, Table 3). The 7.7 and 11.3 /im PAH luminosities are 1.0 ± 0.5 x 10^^ erg s'^ 
and 1.3 ± 0.3 x 10^^ erg s~\ respectively. The 11.3/7.7 PAH fiux ratio is 1.2 ± 0.6, to 
be compared with an average ratio o f 0.3 ± 0.1 for spiral galaxies and AGN host galaxies 



(jSmith et al. 1 120071 : IShi et al. 1 120071 ). The total luminosity in detected PAH features is 



L(PAH) = 1.0 ± 0.1 X 10^2 g-i and the ratio L(H2)/L(PAH)= 0.8 ± 0.1. This ratio is 
much greater than the typical ratio of ~ 1 x 10~^ seen in normal star-forming galaxies where 
both H2 and PAHs are thought to be produced in stellar PDRs ( iRoussel et al. 1120071 ). 



PAH emission excited by UV photons from O and B stars roughly correlates with other 
star- formation rate (S FR) indicators such as Ha luminosity and 24 fim continuum luminosity 
(IRoussel et al. II2OOII ). However, there is a la rge scatter, with PAH luminosity depending on 
metallicity and age of the stellar popu lation (jCalzetti et al. 1120071 ). X-ray radiation from an 
AGN can destroy PAHs (IVoit Ill99ll ). but should have little effect at large distances or in 
regions shielded from the nucleus. With these caveats in mind, we estimate the SFR from 
the luminosity of the 7.7 /im PAH feature to be 7± 3 x 10"^ yr^^ for 3C 326 N. The 24 
Hm luminosity measured with MIPS yields a comparable SFR of 0.14 ± 0.02 Mq yr~^. The 
similarity of these two estimates indicates that the AGN does not dominate the continuum 
at this wavelength (< 50%). The star formation rate is only 2 — 5% of the ~ 3 Mq yr~^ 
SFR of the Milky Way. 

A remarkable aspect of 3C 326 N (and MOHEGs in general) is the lack of starburst 
activity in the presence of a large molecular gas mass. The mean surface density of warm H2 
inside the SL slit is 30 Mq pc~^. In order to suppress star forr aation in a centra l thin disk, the 
total gas surface density must be below the critical density (IKennicutt Ill989l ). An accurate 
estimate of the critical density in 3C 326 N will require a kinematic study of the molecular 
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gas. A large turbulent velocity dispersion or rotational velocity shear would increase the 
critical density for star formation. The molecular gas may have not yet even settled into a 
disk. Efficient star formation and AGN accretion may be suppressed and delayed in 3C 326 
N until the molecular gas has kinematically settled and cooled. 

The PAH and 24 fim luminosities of companion galaxy 3C 326 S yield estimated SFRs 
of 0.20±0.05 Mq yr-i and 0.26±0.03 Mq yr'^ respectively. Thes e SFR estimates are within 



the b road range seen for normal, non-starburst early type galaxies (jCombes. Young, fc Bureau 



20071 ). The SFR in 3C 326 S is significantly larger than the SFR in 3C 326 N. Since no 
emission lines are detected, it is likely that most of the remaining molecular gas is cold and 
conducive to star formation. A Hubble Space Telescope image of 3C 326 S in the F702W 
band shows a dust lane (or nearly edge-on disk) cutting across the co re at PA~ 105°, which 



confirms the presence of a significant dusty ISM (iMartel et al. 1119991 ). The galaxy 3C 326 S 
may be the H2 donor for 3C 326 N, and may have had a much larger molecular gas compo- 
nent which has been largely stripped by tidal interaction. It would be particularly useful to 
measure the masses and surface densities of cold CO (and thereby H2) in both galaxies to 
see if they are consistent with the observed SFRs and to help understand why the warm II2 
emission is so much stronger in 3C 326 N. 



7. Energetics and H2 Heating Mechanisms 

In 3C 326 N, the ratio of mid-IR rotational H2 emission to the 8-70 /im IR continuum 
emission is L{}i2) / L^-jo = 0.17 ± 0.02. The H2 pure rotational lines are therefore a major 
coolant for the warm molecular gas phase of the ISM. The 125 K H2 component has a thermal 
energy content of ~ 4 x 10^^ erg and the cooling time from the observed luminosity in H2 
S(0)-|-S(1) lines is only 9 x 10^ yr. It is therefore necessary to continuously inject energy into 
the molecular gas to maintain the observed temperature. We consider below two possible 
scenarios for heating the H2, i) shock-heating by the inner radio jet and ii) tidally-induced 
infiow or accretion from the companion galaxy 3C 326 S. 



7.1. Radio Jet Heating Scenario 



The radio lobe size and Alf yen speed give an estim ated expansion age of ~ 2 x 10^ yr 



for the Mpc-scale radio source (jWillis fc Strom I Il978l ). The terminal working surfaces of 
the radio jets are in intergalactic space, far outside both galaxies. Any jet-heating of H2 
must be done locally by the inner jet at a radius of < 3.1 kpc. The average jet power is 
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best estimated from the extended radio lobe emission. The minimur n energy content in th e 
radio lobes derived from equipartition arguments is 6.0 x 10^^ erg (jWillis fc Strom Ill978l ). 
Dividing by the radio source lifetime, we estimate a jet kinetic luminosity of 9.5 x 10^^ erg 
s~^. We estimate a smaller kinetic luminosity of 4.8 x 10^^ erg s~^ if the lobe magnetic fields 
have sub-equipartiti on strength and the lobe energy is dominated by relativistic particles 

Jpunsivlbooaboosh . 

The integrated 10 MHz-10 GHz synchrotron luminosity from the radio lobes is 4.7 x 10^^ 
erg s~^, ~ 5% of the estimated jet kinetic luminosity. The synchrotron luminosity of the 3C 
326 N core is much weaker, at 3.5 x 10^^ erg s~^. If we assume that the synchrotron luminosity 
is roughly proportional to the dissipation of jet kinetic energy, then the ratio of core to lobe 
total luminosity indicates that ~ 7 x 10^^ of the jet kinetic luminosity, or ~ 7 x 10^° erg 
s~^ is dissipated in 3C 326 N. (For the sake of illustration, this assumes that relativistic 
beaming is not important in the core. In the likely case that the core emission is beamed, 
this should be divided by the beaming factor.) This is less than one tenth of the luminosity 
required to power the observed H2 emission. Also, jet shocks must be active over a large 
fraction of the total H2 mass in the galaxy to be a viable heating mechanism. However, with 
a large uncertainty in the amount of jet kinetic power dissipated and the unknown spatial 
configuration of H2, this scenario may deserve further consideration. High-spatial resolution 
radio observations are necessary to look for any direct evidence of jet-ISM interaction. 



7.2. Tidally Induced Inflow 



Another, preferred scenario is that 3C 326 N has been disturbed by gravitational in- 
teraction with 3C 326 S. The H2 may have either been stripped from 3C 326 S, or already 
present in 3C 326 N. Tidal forces from 3C 326 S would then induce the observed distor- 
tions in 3C 326 N, in turn giving rise to gravitational torque that drives molecular gas into 
the galaxy center. The gravitational energy of the in-falling gas is then converted into tur- 
bulence and heat via an ac cretion shock or shocks between sub-clumps. The simulations of 
Mihos &: Hernquist I (119961 ) show that the greatest gas inflow rates are achieved for co-planar 
interactions of disk systems ~ 2 x 10*^ yr after the initial encounter due to the growth of a 
central bar over the disk dynamical time scale, but before the final merger. The tidal field 
of the bar drives the inflow. 



T he K-band magn i tudes of 3C 326 N and S are 12.95 and 13.60, respectively (ILilly fc Longair 



19841 : iRawlings et al. I [19901) ■ We estimate stellar bulge masses of 3.1 x 10 Mq and 



1.6 X 1 0^^ Mf7), respectively, us ing an empirical relation between K-band luminosity and bulge 
mass ( iMarconi Sz Hunt II2003I ). The ~2:1 mass ratio qualifies this as a major interaction. The 
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corresponding nuclear supermassive black hole masses (assuming MBn/^buige = 1-5 x 10~^) 
are ~ 5 x 10^ Mq and ~ 2 x 10® Mq, which may eventually result in a major black hole 
merger if the two galaxies merge. The isophotes of 3C 326 N are roughly elliptical at a radius 
of 12'.'2, with PA= 154°. At larger radii (12 — 16"), the isophotes bend to the south to meet 
3C 326 S, perhaps indicating a tidal bridge. The two galaxies have identical radial velocities 
to within the measurement uncertainties (±300 km s~^). 

The 2-body LI Lagrange point is at 24/ cos(^) kpc from the nucleus of 3C 326 N and 
17/ cos(^) kpc from the nucleus of 3C 326 S, where 9 is the inclination of the orbital axis. The 
gravitational potential energy released by 1.1 x lO^M© of H2 falling from LI to the center of 
3C 326 N is 8.4 x 10^^ erg. Dividing by the dynamical (free-fall) time-scale of tdyn = 9.1 x 10^ 
yr, we obtain an H2 inflow rate of 12 Mq yr and an accretion luminosity of 2.9 x 10^^ 
erg s~^. This is sufficient to power the observed H2 emission at a conversion efficiency of 
27%. A massive dark halo will deepen the gravitational potential and allow an even lower 
conversion efficiency. The near-coincidence of the crossing time-scale 2tdyn = 1.8 x 10® yr 
and the age of the radio source suggests the intriguing possibility that the radio jet activity 
is ultimately fueled, at low radiative efficiency, by tidally driven accretion. 

Ground-based adaptive-optics imaging of near-IR H2 rovihrational transitions in 3C 
326 N will probe the spatial distribution of hot H2 and furt her test the scenario of tidall y 



induced accretion. The proposed H2 Explorer mission (H2EX. lBoulanger fc Maillard 1120071 ). 
sensitive to H2 emission over a large field of view, can establish the space density of the general 
MOHEG population and help determine its importance as a stage in galaxy evolution. 



8. Conclusions 

We have discovered extremely luminous molecular hydrogen pure-rotational emission 
lines from radio galaxy 3C 326 N. The continuum emission is relatively weak, so that the H2 
emission contributes 17% of the 8-70 /xm IR luminosity and is a major coolant. We measure 
a warm molecular gas mass of 1.1 X 10^ Mq at temperatures of 125-1000 K within a 3.1 kpc 
radius of the galaxy center. The H2 ortho/para ratios are consistent with thermodynamic 
equilibrium. 

We consider 3 ways to heat the H2 to the observed temperatures and produce the 
extreme H2 emission line luminosities. The LINER AGN in 3C 326 N is most likely not 
luminous enough to power the observed H2 emission via X-ray heating. Jet-shock heating 
could work in principle, but there does not appear to be enough kinetic energy dissipated 
inside the host galaxy by the jet. We favor the idea that accretion shocks from inflow induced 
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by tidal interaction with 3C 326 S may power the remarkable H2 emission spectrum. 

Weak PAH emission features and thermal dust continuum in 3C 326 N and 3C 326 S 
indicate very little star formation activity ~ 0.1—0.3 Mq yr~^. The weak star formation in 3C 
326 N and low accretion luminosity of the AGN are remarkable considering the large amount 
and high surface density of molecular gas. Similar H2 emission seen in some ULIRGs is 
usually accompanied by vigorous starburst or AGN activity. The newly discovered molecular 
hydrogen emission galaxy (MOHEG) phenomenon may represent an important stage in 
galaxy evolution. Star-formation as well as accretion onto the central supermassive black 
hole may be suppressed by high temperatures, turbulence, or velocity shear in the large 
central mass of shocked molecular gas. 

This work is based on observations made with the Spitzer Space Telescope, which is 
operated by the Jet Propulsion Laboratory (JPL), California Institute of Technology (Gal- 
tech) under NASA contract. Support for this research was provided by NASA through an 
award issued by JPL/Galtech. This research has made use of the NASA/IPAG Extragalactic 
Database (NED) which is operated by JPL/Galtech, under contract with NASA. The IRAG 
3-color image was processed using Montage, funded by NASA's Earth Science Technology 
Office, Computation Technologies Project, and maintained by the NASA/IPAG Infrared Sci- 
ence Archive (IRSA). We made use of data from the Sloan Digital Sky Survey funded by the 
Alfred P. Sloan Foundation, Participating Institutions, NSF, NASA, U.S. DOE, Japanese 
Monbukagakusho, Max Planck Society, and the Higher Education Funding Council for Eng- 
land. Thanks to K. Sheth for enlightening discussions on star formation in galaxies. Thanks 
also to T. Jarrett and R.-R. Chary for help with IRAG. 
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Tabic 1. H2 0-0 S(J) Emission Lines 





H2 S(0) 


S(l) 


S(2) 


S(3) 


S(4) 


S(5) 


S(6) 


S(7) 


A(/xm, rest) 


28.22 


17.04 


12.28 


9.66 


8.03 


6.91 


6.11 


5.51 


Flux* 


0.24(0.03) 


0.58(0.05) 


0.39(0.04) 


1.18(0.04) 


0.28(0.06) 


0.86(0.15) 


0.24(0.08) 


0.28(0.08) 




1.9(0.6) 


1.3(0.2) 


0.41(0.08) 


0.80(0.08) 


0.14(0.03) 


0.34(0.10) 


0.074(0.025) 


0.066(0.022) 


Luminosity'^ 


0.48 


1.1 


0.77 


2.34 


0.55 


1.7 


0.48 


0.55 


N/g'^ 


9.2E18 


2.8E17 


3.8E16 


6.9E15 


1.3E15 


4.4E14 


1.5E14 


2.7E13 



"^Flux (10"^'' erg cm~^). 
^Equivalent Width (/^m, rest). 

Luminosity (lO"'^ erg s^^), at a luminosity distance of 407 Mpc. 
^Column Density (cm~^) of upper level divided by its statistical weight, assuming a 3'.'7 x 3'.'7 source. 
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Table 2. Forbidden Emission Lines 





[Fe II] 


[SIV] 


[Ne II] 


[S III] 


[OIV] 


A(/im, rest) 


5.34 


10.51 


12.81 


18.71 


25.89 


Flux'^ 


0.20(0.08) 


0.08(0.04) 


0.40(0.04) 


0.14(0.07) 


0.25(0.03) 


Flux/[Ne II] 


0.5(0.2) 


0.2(0.1) 


1.0 


0.4(0.2) 


0.6(0.1) 



"^Flux (10-1^ erg s'^ cm-^). 
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Table 3. Aromatic Emission Features (PAHs) 



A(/im) 


6.7 


7.7 


10.7 


11.3 


13.6 


15.9 


17 


18.9 


3C 326 N 


1.13(0.48) 


0.53(0.26) 


0.22(.09) 


0.64(0.13) 


0.98(0.18) 


0.21(0.11) 


1.15(0.27) 


0.33(0.18) 


3C 326 S 




1.61(0.42) 




1.19(0.20) 











"^Flux (10-1^ erg s'^ cm-^). 
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Table 4. MIPS 24, 70, and 160 /xm Photometry 





F^(24 /im)^ 


F,(70 nm}^ 


F^(160 /xm)'^ 


L(24 /im)'' 


L(70 /im) 


L(70)/L(24) 


3C 326 N 


0.52±0.07 


6.12±0.65 


< 34 


1.3±0.2 


5.2±0.6 


4.0±0.8 


3C 326 S 


1.05±0.12 


4.33±0.60 


< 34 


2.6±0.3 


3.7±0.5 


1.4±0.3 


bridge'' 




4.40±0.60 






3.8±0.5 





^Fhix density (mjy) in the 24 /xm band, measured in a 4'/9 radius aperture and multiplied by 

an aperture eorreetion of 2.0. 

•^Flux density (mJy) in the 70 iim band, measured in a 9'.'84 radius aperture and multiplied by 

an aperture correction of 1.5. 

■^The 3(7 upper limit to the 160 /im flux density (mJy), measured in a 20" radius aperture. 

Luminosity vL^{\{)'^^ erg s^^), at a luminosity distance of 407 Mpc. 

''This source (only in the MIPS 70 /im image) is found at roughly the location where the tidal 
bridge appears to connect to 3C 326 N (Fig. 5). 
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Table 5. H2 Model Parameters 



T(K) 


M{H2T 


Ortho/Para 


^crit 


./max 


125 


1.1E9 


2.1 


5.3E2 


1.5E-3 


400 


1.1E7 


3.0 


1.1E4 


7.0E-7 


1000 


6.4E5 


3.0 


4.4E5 


l.OE-9 



'*H2 mass {Mq). 
''Critical density (cni~^). 

'^Maximum filling factor under LTE conditions. 
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Table 6. Strong H2 Emitters 





M(cold Hs)"* 


M(warni 


L{H2r 


L{H2)/L{IRY 


Refs. 


Zw 3146 BCG 


~ 1 X 10" 


~ 1 X 10^° 


2.1 X 10*3 


0.014 


1 


NGC 6240 


1.1 X 10^° 


1.6 X 10^ 


4.6 X 10^2 


0.003 


2,3 


3C 326 N 




1.1 X 109 


8.0 X 10^1 


0.04-0.1 




Antennae 


9.6 X 10^ 


4.9 X 10^ 


1.7 X lO^^i 


0.0009 


4,5,6 


Stephan's Quintet 




> 3.4 X 10^ 


> 8.4 X 10*° 


0.34 


7 



'^Cold H2 mass derived from CO measurements. 
''Warm H2 mass measured from Spitzer or ISO spectra. 
■^112 Luminosity (erg s~^). 

'^The IR luminosity is integrated over 8-1000 /xm. 



References. — (1) Egami et al. 2006; (2) Lutz et al. 2003; (3) Armus et al. 2006; (4) Haas 
et al. 2005; (5) Young et al. 1995; (6) Klaas et al. 1997; (7) Appleton et al. 2006. 
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Fig. 1.— Spitzer IRS SLl 8.5-13.4 /im (rest) 2D spectra of 3C 326 N (top) and 3C 326 S 
{bottom). The bright H2 S(3) 9.66 /im emission line in 3C 326 N is spatially and spectrally 
unresolved. The H2 S(2) and Ne II lines are also seen. The southern galaxy was serendip- 
itously placed precisely on the SL slit (PA= 178°). Its 11.3 /xm PAH emission feature is 
indicated. 
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0-0 S(7) S(6) S(5) S(4) S(3) S(2) S(1) S(0) 

I I I I I I I 3C 326 N 
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45678910 20 30 
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Fig. 2. — Top: Spitzer IRS spectrum of 3C 326 N. Dashed lines-SLS and LL3 bonus orders. 
Extremely strong H2 emission lines of high equivalent width dominate the spectrum. Weak 
PAH features indicate a very low level of star formation in spite of the large amount of H2. 
Forbidden ionic emission line ratios are consistent with a weak LINER AGN. Bottom: PAH 
emission is detected in 3C 326 S, but no H2 or ionic forbidden hues. 
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Fig. 3. — Spitzer IRAC 3.6, 4.5, 5.8, and 8.0 fim images of 3C 326. North is at top. East to 
the left. The first two bands are dominated by stellar emission, while the second two have 
contributions from dust and H2. The galaxies 3C 326 N and 3C 326 S are the brightest 
two sources in these images. A third source at PA = 268° and a fourth 4'.'3 away from 3C 
326 N at PA = 111° (seen protruding from the left of the galaxy core in the 4.5 //m image) 
appear to be unrelated galaxies seen in projection (see §3.3). Note the apparent tidal bridge 
between 3C 326 N and S and a possible tidal tail to the NW in the 3.6 and 4.5 /im images, 
evidence of an ongoing galaxy collision. The galaxy cores are barely resolved at 5.8 /xm, and 
3C 326 S is nearly as bright as 3C 326 N at 8.0 /xm. 
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Fig. 4. — Spitzer IRAC 3-color image of 3C 326. Blue = 3.6 /im, green=4.5 /xm, red=5.8 
/iin. Foreground stars and starlight-dominated galaxies appear blue-green, while dust or 
H2 emission appears red. The brightest galaxy and AGN host 3C 326 N has an irregular 
morphology because of its interaction with 3C 326 S. Note the possible tidal tail extending 
to the NW. 
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Fig. 5. — Left: Spitzer MIPS 24 /xm flux contours overlaid on the IRAC 3.6 /xm image of 3C 
326. North is at top, East to the left. The contour levels are 0.3, 0.4, 0.5, and 0.6 MJy/sr 
Right: MIPS 70 iim contours overlaid on the 3.6 iim image. The contour levels are 0.45, 
0.65, 0.75 and 0.85 MJy/sr. 3C 326 S is brighter at 24 //m, while 3C 326 N is brighter at 70 
fim. A third 70 fim source appears between the two galaxies, but not in the IRAC image. 
It may be cool dust emission associated with the tidal bridge, or an unrelated background 
source. 
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Fig. 6. — Spectral energy distributions of 3C 326 N {open triangles and solid line), 3C 
326 S (squares), and 3C 326 radio lobes (diamonds). MIPS 24 and 70 /im points indicate 



warm and cold dust components. MIPS 160 fim and ground-based 870 fim (iQuillen et al. 



20031 ) 3cr upper limits (inverted tri angles) apply to both g alaxies individually. J, H, and 



K-band near-IR photometry are bv iLillv &: Longair I (19841). Th e radio core fluxes and the 



K-band photometry of 3C 326 N are from iRawlings et al. I ( 



19901 ). Radio lobe data are from 



WiUis fc Strom I (119781 ) and the NASA/IPAC Extragalactic Database (NED). The dashed 
line represents the Rayleigh- Jeans tail of starlight connecting the K-band photometry to the 
IRS spectrum of 3C 326 N. The dotted lines are 25 and 58 K blackbodies, representing upper 
and lower limits to the 70-1000 fim emission. 
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Fig. 7. — H2 excitation diagram for 3C 326 N. Upper-level column densities normalized 
by their statistical weights are plotted as a function of upper-level energy {E^/k-o) for each 
pure-rotational transition. The best-fitting model {solid line) has 3 temperature components 
{dashed lines): 125, 400, and 1000 K. The ortho/para ratios are set at the respective thermal 
equilibrium values (Table 5) . Note that the y-intercept for each model component gives the 
average H2 ground-state column density integrated over the SL slit. 



